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Abstract. The cloud of cold atoms obtained from a magneto-optical trap is known to exhibit two types of 
instabilities in the regime of high atomic densities: stochastic instabilities and deterministic instabilities. 
In the present paper, the experimentally observed stochastic dynamics is described extensively. It is shown 
that it exists a variety of dynamical behaviors, which differ by the frequency components appearing in the 
dynamics. Indeed, some instabilities exhibit only low frequency components, while in other cases, a second 
time scale, corresponding to a higher frequency, appears in the motion of the center of mass of the cloud. 
A one-dimensional stochastic model taking into account the shadow effect is shown to be able to reproduce 
the experimental behavior, linking the existence of instabilities to folded stationary solutions where noise 
response is enhanced. The different types of regimes are explained by the existence of a relaxation frequency, 
which in some conditions is excited by noise. 



PACS. 32.80.Pj Optical cooling of atoms; trapping 
nonlinear dynamical systems 



1 Introduction 

Magneto-optical traps (MOT) produce clouds of cooled 
atoms at temperature as low as the /iK. The achievement 
of such clouds opened many perspectives, not only in the 
field of fundamental atomic physics, as e.g. in the domain 
of the atomic dynamics or the quantum chaospQ, but also 
leads to several potential applications, in particular the 
improvement by several orders of magnitude of the atomic 
clockjjj. The MOT is also the first stage to produce lower 
temperatures, in particular to obtain Bose-Einstein con- 
densates ■ Although the use of MOTs is relatively well 
mastered, some details of the experimental setup remain 
empirical, because of the existence in the cloud of insta- 
bilities that are not well understood. Some studies showed 
that when the trapping beams are misaligned, the cloud 
may be spatially altered and become unstable. In partic- 
ular, ring-shaped clouds and chaos have been observed, 
and attributed to a vertex force However, the main 

instabilities encountered in the experiments concern well- 
aligned MOTs. In that case, the cloud, which is usually 
more or less ellipsoidal, has a complex irregular shape, 
with an inhomogeneous atomic density (figGJ. This shape 
may not be stable, and changes as a function of time: a 
typical example of these instabilities is shown in fig. ^ 

There is a deep interest to identify the nature of these 
instabilities, in the aim to control them, and possibly to 
take advantage of them. To illustrate these points, let us 
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Fig. 1. Sequence of snapshots showing the time evolution 
of the unstable atomic cloud. Snapshots are presented in the 
chronological order, each one being separated by 120 ms. 



remember that instabilities may originate in many mech- 
anisms, which can be classified in two families: stochastic 
or deterministic. In the latter, called also deterministic 
chaos, the dynamics is described by a set of deterministic 
equations. When the number of equations is small (low- 
dimensional deterministic chaos), it has been shown that 
the dynamics can be controlled to reach various states 
that are not accessible otherwise, as unstable states|H] or 
periodic behaviors^- It is also well known that the de- 
scription of such a dynamics gives access to numerous pa- 
rameters, which are not measurable when the behavior 
is stationary 8 . The cost for these informations is that 
there is no other way to reduce the instabilities. On the 
contrary, if the origin of the instabilities is stochastic, i.e. 
due to external noise, a reduction of instabilities is ob- 
tained by reducing the noise, but the instabilities have no 
physical meaning, in that sense that they are not intrinsic 
to the physics of the system. 
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In 0, a first type of instabilities, arising for low beam 
intensities, has been depicted. It was showed that these 
so-called "stochastic instabilities" are induced by the ab- 
sorption of the trapping beams in the cloud. A model tak- 
ing into account this so-called shadow effect showed that 
from a dynamical point of view, instabilities arise through 
a stochastic resonancelike phenomenon, namely the coher- 
ent resonance, linked to a Hopf bifurcation in the station- 
ary solutions of the MOT. It is well known that such a 
bifurcation usually leads to periodic instabilities, and in- 
deed, a recent study evidenced such instabilities, purely 
deterministic [IQj . However, a correct description of these 
self-oscillations required to modify the model in j2j , in par- 
ticular to take into account the spatial distribution of the 
cloud. 

In the present paper, the results presented in [5] are 
detailed, and extended to a larger range of parameters, 
where new regimes appear, in particular a behavior with- 
out resonance frequency, contrary to the instabilities de- 
scribed in [Hj . Then the modified model introduced in JOj 
is accurately described, and the relative domains of ap- 
pearance of stochastic and deterministic instabilities are 
discussed. The difference appears to be the result of the 
preeminent role of noise in some parameter range, in par- 
ticular in the vicinity of the bifurcations. The stochastic 
instabilities predicted by the model are studied, and they 
are shown to be in agreement with those experimentally 
observed. In particular, an interpretation of the two types 
of stochastic instabilities, with and without resonance fre- 
quency, is found. 

The paper is organized as follows. After this introduc- 
tion, the section|21describes the experimental setup. Then, 
a detailed analysis of the experimental observations of the 
stochastic instabilities is presented (section showing in 
particular the existence of two types of instabilities. Sec- 
tion|Uis devoted to the construction of a simple ID model, 
already presented in |10j in a less detailed way. In section 
the stationary solutions of the model are discussed. Fi- 
nally, in section^ the noise induced dynamics instabilities 
is studied and compared to the experimental observations. 



2 Experimental set-up 

We work with a Cesium-atom MOT in the usual <7+ — cr_ 
configuration, with three arms of two counter-propagating 
beams obtained from the same laser diode. The waist 
wt of the trap beams may be varied from typically 3 to 
10 mm. Two configurations are possible: in the first one, 
all six beams are independent, by opposition to the sec- 
ond configuration, where counter-propagating beams re- 
sult from the reflection of the three forward beams. In 
the last case, the intensity asymmetry resulting from the 
absorption of the forward beam by the cloud, generates 
a center-of-mass motion, while in the first case, instabili- 
ties are characterized by symmetrical bursts on the cloud 
shape, much more difficult to measure. However, as the 
nonlinearities involved in both cases are the same, we ex- 
pect that the dynamics will be fundamentally of the same 



nature, and thus we choose the configuration with retro- 
reflected beams. 

A full description of the unstable dynamics of the atomic 
cloud will be presented in the next section. However, to 
make easier the understanding of this paragraph, let us 
depict them briefly. As shown in the introduction (fig.^l, 
instabilities consist in a deformation of the spatial atomic 
distribution, leading to fluctuations of the shape of the 
cloud. Therefore, the relevant dynamical variables allow- 
ing us to describe instabilities, could be the shape of the 
cloud (i.e. for example the local velocities and atomic den- 
sities in the cloud) . This type of description corresponds to 
a high dimensional model, associated with partial differen- 
tial equations. Here, for the sake of simplicity, we choose 
to limit our description to the center of mass (CM) lo- 
cation r, and the total number of atoms n in the atomic 
cloud. This allows us to model the system with ordinary 
differential equations, and reduces the dimension to seven, 
and even three in a ID model. As it is shown in the fol- 
lowing, the use of this description appears to be sufficient 
to understand the main mechanisms of the instabilities. 

To measure r and n, we used two 4-quadrant photodi- 
odes (4QP) forming an orthogonal dihedral and measuring 
the fluorescence of the cloud. The differential signal of the 
4QP allows us to monitor the motion of the CM. Using 
only two 4QPs does not allow us to reconstruct the actual 
3D motion, but gives access to projections of this motion 
on two different axes. This prevents the measure from line- 
of-sight effects due to the optical thickness of the cloud. 
We checked that whatever the type of dynamical behav- 
ior, the motion components r recorded by both 4QP have 
the same properties and are qualitatively identical. In the 
following, the "CM motion" refers to a component r of 
this motion recorded by one of the 4QPs. The second dy- 
namical variable, namely the number of atoms inside the 
cloud, is deduced from the total signal received from the 
4QPs. In addition to the 4QPs, two video cameras monitor 
the shape of the cloud. Because of their poor resolution as 
compared with the 4QPs, the video cameras are not used 
to record the dynamical variables. They have been essen- 
tially used in the first stages of the experiment, to control 
that there is no discrepancy between the shape dynamics 
and the CM dynamics. 

Instabilities depend on the MOT parameters. Among 
these parameters, some are easily controllable, and will be 
referred in the following as control parameters. They are 
the detuning Aq of the MOT, the magnetic field gradient 
G, the MOT beam intensities I\ and the repumper laser 
intensity I re p- Other parameters cannot be considered as 
control parameters, because they are not easily control- 
lable or measurable in our experimental setup. Among 
these parameters, the alignment of the MOT beams ap- 
peared to be crucial in the experiment. We limit the present 
study to the case where beams are aligned. When mis- 
aligned, most of the dynamical characteristics of the cloud 
change qualitatively |5] . The MOT beam waists play also a 
main role on the dynamics, but in practice, they cannot be 
changed easily independently from the other parameters, 
in particular I\. However, two different values have been 
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used in the experiments, and their effects on the dynam- 
ics are well understood. Finally, the vapor pressure in the 
cell has probably also a large influence on the dynamics. 
Unfortunately this parameter is not easily measurable in 
our experimental setup. Moreover, as it is shown later, it 
does not appear explicitly in the model, but its impact on 
the dynamics may be estimated through the equilibrium 
population of the cloud. The parameter ranges explored 
in the present experiment are summarized in Tab. ^ 

Table 1. Range of the parameters used in the present exper- 
iment. G is the magnetic field gradient, I\ is the intensity of 
the forward beam, Ao is the detuning and wt is the trap laser 
beam waist. 7 S is the saturation intensity ( I s = 1.1 mW) and 
F is the natural width of the transition. The last column indi- 
cates the default parameter values used to obtain the results 
reported in the present paper. 

range default set 
G (Gem" 1 ) G < 14 14 
7+ = I /I s 4 < 7+ < 20 6 
WT{mm) 3 < wt < 10 3 
A A < -0.5 



3 Experimental results 

Instabilities have been described in [H] for a given set of pa- 
rameters. In the following, we extend this description for 
the whole range of parameters where stochastic instabili- 
ties appear. A first fundamental control parameter is the 
MOT beam intensity, whose value determines the type of 
observed instabilities. At low MOT beam intensities, typ- 
ically less than 1075 (Is = 1.1 mW/cm 2 is the saturation 
intensity), the cloud exhibits S instabilities (S stands for 
Stochastic). For MOT beam intensities larger than 1075, 
C (for Cyclic) instabilities appear JOj • As the aim of the 
present paper is to discuss about the stochastic instabili- 
ties, we keep a detailed presentation of the deterministic 
C instabilities for another paper. However, the two types 
of instabilities cannot be completely separated, and the 
domains of appearance of both types of instabilities will 
be discussed. 

As discussed above, S instabilities appear at low MOT 
beam intensities, and are characterized by large fluctua- 
tions of the shape of the cloud appearing in a limited range 
of the parameters. From the experimental point of view, 
this last point is essential, because it is at the origin of the 
introduction of the concept of instabilities of the MOT. 
Indeed, if the noisy dynamics is the same whatever the 
parameter values, it is clear that the problem becomes as 
trivial as the reduction of technical noise in an experiment. 
On the contrary, the motion of the MOT grows suddenly 
in a narrow range of the parameters, as it is illustrated 
in fig. |5J where the amplitude Ar of the r fluctuations 
is represented as a function of the control parameter Aq. 
Far from the unstable area, e.g. at large \Ao\, the cloud is 
stable in shape and density. When the control parameter 



is tuned to the unstable area, instabilities appear progres- 
sively, and the amplitude grows until a maximum Z\r max 
at Ziomax — —1-95 (Aq is expressed in units of the natu- 
ral width r of the atomic transition). There is no abrupt 
boundary between the stable and unstable areas: the lim- 
its given below corresponds to Ar — Z\r max /10. In this 
case, for the set of parameters of fig. [3 instabilities appear 
for —2.5 < Aq < —1.7 (for Aq > —1.7, the population in 
the cloud vanishes). The unstable range depends on the 
other parameters, such as the beam intensities, the vapor 
pressure or the magnetic field gradient. For example, for a 
less populated cloud, due e.g. to a different vapor pressure 
in the cell, all other parameters being the same as previ- 
ously, instabilities will appear at smaller detuning. How- 
ever, whatever the parameters used in the experiments, in 
the range given in Tab. ^ we have — 3 < Z\o m ax ^ — lj 
and instabilities never occur on a range larger than 1. 




-3.0 -2.8 -2.6 -2.4 -2.2 -2.0 -1.8 -1.6 



A 

Fig. 2. Evolution of the signal amplitude Ar as a function of 
the detuning Ao. Parameters are those of the default set of 
Tab. [I] 

Fig. |2| shows a typical unstable behavior of r. It ap- 
pears as an erratic signal, with a flat spectrum (fig. 
Low frequency significant components appear typically for 
values smaller than v n w 2 Hz (fig. , and the dynam- 
ics is essentially along the first bisector of the three for- 
ward beams. The behavior of r and n are similar, with a 
cross correlation coefficient larger than 0.8. To determine 
if these instabilities have a deterministic origin, several 
tools are offered through the nonlinear dynamical analy- 
sis of the time series. As the MOT is dissipative, a de- 
terministic dynamics should have an attractor, which can 
be reconstructed easily for a low dimensional dynamics. 
The result, not presented here, appears as a set of ran- 
domly distributed points: in particular, it does not present 
any fine structure. Poincare section and ID maps confirm 
this absence of order in the dynamics. This could be due 
to a lack of resolution of the measures, but the general 
shape of the trajectories rather suggests that the behav- 
ior is stochastic (or chaotic with a high dimensional dy- 
namics). Because this behavior appears to be a stochastic 
dynamics with only low frequency components, it will be 
referred in the following as Sl instabilities. 
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Fig. 3. Experimental record of the time evolution of a compo- 
nent of the CM location of the atomic cloud. Experimental pa- 
rameters are the default ones given in Tabled with Ao = — 2. 
The mean cloud population is 1.5 x 10 8 atoms and the cloud 
size is 1 mm. 
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Fig. 4. Power spectra corresponding to the behavior illustrated 
in fig. El In (a) CM location, and in (b) cloud population. The 
scales are linear. 



In some situations corresponding to given ranges of pa- 
rameters .5], the r dynamics is altered by the appearance 
of spontaneous large amplitude oscillation-like bursts. As 
illustrated in fig. the signal inside the bursts is not pe- 
riodic, although it is clearly dominated by a given fre- 
quency. These bursts are relatively scarce, and the global 
shape of the signal remains that of Fig. |3l The bursts 
appear in fact as the most spectacular consequence of a 
deeper change of the dynamics, which is the appearance 
of a second characteristic time. This appears clearly in 
the spectrum of r (Fig. HJi) as a peak centered at a fre- 
quency v r , which depends on the parameter values, but 
ranges typically between 10 and 100 Hz. The other char- 
acteristics of the dynamics are not modified by the exis- 
tence of bursts. In particular, the spectrum still exhibits 
the low frequency component below v n , and the bursts do 
not appear on the n behavior (Fig. andHjJj). The low 
frequency component dynamics of r and the dynamics of 
n remain correlated. This behavior will be referred in the 
following as Sh instabilities. 

When the trap beam intensity I\ is increased, 5" insta- 
bilities still exist, but they are progressively superseded by 
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Fig. 5. Experimental record of the time evolution of a com- 
ponent of the CM location of the atomic cloud in the case of 
the S behavior. Experimental parameters correspond to the 
default set given in Table0 with Ao = —1.5. The mean cloud 
population is 1.5 x 10 8 atoms and the cloud size is 1 mm. 
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Fig. 6. Power spectra corresponding to the behavior illustrated 
in fig. In (a) CM location, and in (b) cloud population. The 
scales are linear. 



C instabilities. These instabilities differ drastically from S 
ones^J: they can be either periodic or erratic, but in both 
cases, they are cyclic, and the motion amplitude is much 
larger. However, as for S instabilities, C instabilities exist 
in a limited range of Aq, which value depends on the other 
parameters. But it is systematically between Aq = —3 and 
resonance, and the maximum detuning range is of the or- 
der of 1. 
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As I\ is increased, the disappearance of S instabilities 
occurs progressively, in favor of C instabilities. For inter- 
mediate values of I\, both types of instabilities exist. Their 
typical distribution versus Aq is illustrated in Fig. far 
from resonance, the cloud is stable; as the resonance is ap- 
proached, S instabilities appear for a detuning Aq = A\. 
Then C instabilities appear in A2 > A\. If the detuning is 
still increased, C instabilities disappear in A3 at the ben- 
efit of a stable behavior. Finally, the cloud vanishes in A4. 
As I\ is increased, the width 8\2 = A2 — A\ decreases in 
favor of the interval ($23 = A3 — A2, while the total unsta- 
ble interval 613, = A3 — A\ remains more or less constant. 
When C instabilities merge for I\ = AI$ , they appear on a 
narrow interval 823 ^ 0. This interval increases rapidly un- 
til I\ = 7.5/5 and £23 = 0.8. For 1\ > 7. 5/5, S23 increases 
more slowly, to reach the value of 623 — 1 in I\ = 201$, 
where S instabilities have completely disappeared. 
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Fig. 7. This figure illustrates the evolution of the behavior as a 
function of the detuning for I\ = 6.8 and I rep = 1.5 mW/cm 2 . 
The full line reports the population, while the dashed lines 
separate the domain of different behaviors: st. stands for stable, 
S for S instabilities and C for C instabilities. 



To conclude this section, let us summarize the main 
properties of the instabilities: they appear to be stochas- 
tic at a time scale larger than 0.5 s. They may exhibit a 
second dynamical component with time scale smaller than 
10~ 2 s and acting only on r. The low frequency component 
corresponds to a CM motion along the first bisector of the 
three forward beams. Finally, they appear only in a lim- 
ited range of the parameters, suggesting a nonlinear origin 
of the instabilities. This last point was also suggested by 
different tests showing that beam phase fluctuations are 
not at the origin of instabilities 9 . 



4 Model 

To understand the dynamics observed in the experiments, 
we build a phenomenological model taking into account 
the shadow effect. The aim here is not to model as finely 
as possible the experimental system, but on the contrary 
to make a model as simple as possible, enlightening the 
fundamental mechanisms leading to the instabilities. 

The model that we used here has already been de- 
scribed in JOj. It is based on the shadow effect induced 
by the intensity gradients produced by the absorption of 



the trapping laser beams in the cloud |Hlllj . The shadow 
effect generates a force compressing the cloud. In the case 
where the counter-propagating beams result from the re- 
flection of the three forward beams, this compression is 
accompanied by a displacement of the CM along the first 
bisector of the three forward beams: indeed, the backward 
beams are less intense than the forward ones because of 
the absorption in the cloud, and so the latter literally push 
the cloud off its equilibrium location. The main role of the 
shadow effect in the instabilities is suggested in the exper- 
iments (i) by the observed correlation between the n and 
slow CM dynamics and (ii) by the fact that the slow dy- 
namical component of r displaces the cloud along the first 
bisector of the three forward beams. The fundamental role 
of the shadow effect in coupling the population and the 
CM position of cold atom clouds was already reported in 
|12| . In that experiment, bistability in the population dy- 
namics was observed by perturbing the atomic cloud with 
a highly-focalized laser beam, and shadow effect is shown 
to be the dominant nonlinearity. 

We built a ID model taking into account the shadow 
effect in a MOT where counter-propagating beams result 
from the reflection of the three forward beams. As we are 
interested here in the collective motion, we use as variables 
the location z of the center of mass of the cloud along the 
unique axis z of the system, and the number of atoms n 
inside the cloud. The origin of z coincides with the "trap 
center", that is, the zero of the magnetic field. Thus, the 
motion of z may be described by the equation: 



dt 2 



1 

M 



(i) 



where M is the mass of the cloud and Ft the global 
force exerted on the atoms by the two counterpropagat- 
ing beams. To evaluate Ft, we assume a multiple scatter- 
ing regime, i.e. a constant atomic density p in the cloud. 
The atoms are distributed between z m j n and z max , and the 
quantity Az = z max — z m j n is the longitudinal size of the 
cloud. We introduce a constant cross section S, allowing 
us to connect Az to n: 



n = pSAz 



(2) 



The repulsive force induced by the multiple scattering is 
an internal force, and thus does not affect directly the 
center of mass motion. It is also the case for the attrac- 
tive force induced by the shadow effect, which acts on 
Ft only through the asymmetry between the forward and 
backward beams. The forward beam is polarized er + and 
has an input intensity I\. After the crossing of the cloud, 
its output intensity is I2. This is also the input intensity 
of the backward beam, which is polarized ct_. After the 
crossing of the atomic cloud, the output intensity of the 
backward beam is I3. We have obviously I\ > I2 > h- Ft 
is proportional to the number of photons absorbed per 
unit time, i.e. S(I\ — Izj/huiL for the forward beam and 
S(l2 — h)/f luJ L for the backward beam (Hujl is the energy 
of one photon, cul being the laser angular frequency). Ft 
is obtained by multiplying the difference between both 
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F T = -{h-2h+h) 
c 



(3) 



To evaluate I2 and I3, we need to solve the propagation 
equation of light inside the cloud. To simplify the calcula- 
tions, we choose to model a F g = — > F e = 1 transition: 
this is the simplest one leading to a magneto-Doppler ef- 
fect. This choice does not allow us to reproduce the sub- 
Doppler effects, neither the photon redistribution between 
the two waves, but as shown in |10| . it allows us to explain 
most of the experimental behaviors. 

The z axis is chosen as quantification axis, and we 
introduce the ground state \g) and the three Zeeman sub- 
levels of the excited state |eo) et \e±). Because of the ki- 
netic momentum conservation, |eo) is not coupled to the 
light field, and thus, the model is reduced to a three level 
system, {\g), |e_), |e + }}. The interaction with light is gov- 
erned by the Rabi frequency f2± and the effective detun- 
ings A±, taking into account the Doppler et Zeeman vari- 
ations: 



A, 



1 



(AoTkvTuj'Bz) = -(A ±2S) 



(4) 



where k is the wave vector of light and w' B the Zeeman 
shift, measured in angular frequency by unit of length. 
The value of the atomic variables is given by the station- 
ary solution of the optical Bloch equations for the density 
matrix a: 

cr = -r [H , a } + & re [ ax (5) 
in 

The hamiltonian H of coupling with light may be written 
as a function of A± and Q± in the {\g), |e_), |e+)} basis: 



Q-/2 Q+/2 
II = h I Q*_/2 -A- 

ay 2 -a, 



(6) 



The relaxation rates a re i ax are different for the excited 
states, the optical coherences and the ground state: 



{&ee)relax — F(7 ee 
(&eg)relax — ^^eg 
(<7gg)relax = ~ r((T ++ + <T ) 



(7a) 
(7b) 
(7c) 



The time scale of the internal variables is -T -1 , and 
so is much shorter than that of the external variables z 
and v. Thus the stationary solution of Eq. J^J can be 
used. The solution for the populations of the excited states 
II± = u±± can be obtained analytically. On the other 
hand, as the redistribution of photons between the beams 
is forbidden, II± is directly proportional to the absorp- 
tion rate of the beams, and the intensity gradient can be 
written: 

dl± 

dz 



(8) 



By injecting the analytical solution of II± in this equa- 
tion, we obtain the following equations of propagation: 



dl 4 

dz 

dz 



-Ffrwp 



— — = rhiop 



1 - A.LUa 



1 - Aa.U A 



where: 



I_± 

'2 7 



[a + p + + p- ) + 



0± = 



+ A + A_ + — {jjl - 1) 
I ± (A ± 6) (~ + A (p+ - ^ 



A{\ + 5 2 ) 
H = p+ + (i- 

6 = f3 + (3_ + a (f3+p + + (3-[i-) 
A± = 2(3 T + a(2p ± -n T ) 
C± = -0± ±a(fj,+ - //_) 



(9a) 
(9b) 

(10a) 

(10b) 

(10c) 
(lOd) 

(lOe) 

(lOf) 
(10g) 
(10b) 
(lOi) 



These equations appear as the ratio of two polynomials 
of high order, and an intuitive interpretation of this re- 
sult is difficult. However, a numerical resolution of these 
equations leads to a rigorous evaluation of Ft- 

The cloud population dynamics is modeled by a "feed- 
loss" rate equation [T3] : 



dn 
~dt 



= B(n e — n) 



(11) 



where we have introduce the population relaxation B and 
the atom number in the cloud at equilibrium n e , which is 
linked to the loading rate L by L — Bn e . n e is assumed 
to depend on the CM location, to take into account the 
losses variation when the cloud moves from the trap cen- 
ter. We do not know the exact form of this dependence, 
because this variation may have several origins. However, 
we can suppose that the main contribution comes from 
the transverse distribution of the trap laser beam, which 
is gaussian. For sake of simplicity, we keep only the first 
terms from its Taylor's series, i.e. the quadratic term. In 
fact, we checked that the exact dependence of n e does not 
change drastically the results given below. So we write: 



n 



1-1^ 

zo 



(12) 



where no is the equilibrium cloud population at the trap 
center and zq a characteristic length. Atoms in \z\ > \zo\ 
are considered as lost: this is taken into consideration 
when n is deduced from Eq. [21 

Finally, we introduce the reduced variables Z = z/zq, 
V = v/v r and N — n/no, where v r is the recoil velocity 



Daniel Hennequin: Stochastic dynamics of the magneto-optical trap 



7 



(v r = hk/m), and we obtain the following autonomous 
system of equations: 



dZ 
~dt 
dV 
~dt 
dN 
~dt 



V— 

zo 

1 



= B(1-Z' i -N) 



(13a) 
(13b) 
(13c) 



Most of the theoretical parameters arc the exact coun- 
terpart of the experimental parameters, as e.g. the mag- 
netic field gradient or the beam intensities. In this case, 
we used in the model the same values as those of Tabled 
It is not the case for all parameters, either because of the 
simplicity of the model or because they cannot be mea- 
sured easily in the experiment. For example, no has not a 
simple experimental counterpart, but depends on several 
experimental parameters, as e.g. the repumping laser in- 
tensity or the vapor pressure in the cell. For this reason, 
a large interval of no values has been used for the simu- 
lations (Table |2J) . The density p and the cross sectional 
area S, which play the same role, have a meaning only 
in the context of a ID model, while they correspond to 
variables in the experiments. They are fixed in the sim- 
ulations at experimental averaged values, and they have 
been varied on a wide range to check that their value is 
not critical. Finally, the parameter z$ has not exact ex- 
perimental counterpart, as it is linked to both the trap 
beam waist and intensities. Indeed, the relevant size for 
the atoms is not the beam waist, corresponding to an in- 
tensity decreased by a ratio e -2 as compared to the center 
of the beam, but rather the location where the local beam 
intensity decreases under Is- This value is much larger 
than wq for intense beams. 



Table 2. Parameters used in the numerical simulations. The 
range corresponds to the interval explored numerically, while 
the different sets refer to most of the results presented in this 
paper. 





range 


set #1 


set #2 


G (Gem" 1 ) 


14 


14 


14 


B (s- 1 ) 


3 < B < 30 


3 


3 


h 


2 < h < 30 


25 


10 


p (cm -3 ) 


10 10 < p < 3 x 10 10 


2 x 10 10 


2 x 10 10 


S(m 2 ) 


10" 6 < S < 3 x 10 -6 


10" 6 


10" 6 


z (m) 


10" 2 < z < 3 x 10" 1 


3 x 10" 2 


3 x 10~ 2 


no 


10 7 < n < 10 9 


10 8 


2 x 10 7 


A 


5 < A < 


-1.5 


-0.23 



In order to check if the model can be more simplified, 
we tried to reproduce the unstable dynamics of the cloud 
with the first terms of the Taylor's series of Eqs. I|9a|l et 
H9b(l . but we needed to keep several orders and did not ob- 
tain simpler equations. Another possible approximation 
concerns the different terms appearing in A± , in Eq. ^ 
Using the values given in table II, one sees easily that the 
Aq term and the Zeeman shift are of the same order of 



magnitude, while the Doppler shift goes to zero at equi- 
librium, but can be the largest term out of equilibrium. In 
these conditions, no approximation is possible. Therefore, 
we use the equations (II 3D in the calculations. 

To perform the comparison between the experiments 
and the present model, we need to study the behavior of 
the system when noise is added. Noise can be added in Eqs 
1131 on any parameter, and as we were not able in our ex- 
periments to identify clearly the main source of noise, we 
tried theoretically several parameters, as the beam inten- 
sity Ii or the equilibrium population uq. In this case, the 
stochastic model corresponds to Easll3lwhere the parame- 
ter Ii (resp. no) is replaced by (1 + C) Ii (resp. (1 + () n ), 
where £(t) is the noise component. We also simulated sev- 
eral types of noise: gaussian white noise, but also colored 
noise with different distributions of the frequencies. All 
configurations give identical results: on the one hand, the 
choice of the noisy parameter is not critical; on the other 
hand, the spectrum of noise does not alter the response 
spectrum, except obviously that it depends on the rela- 
tive weight of the different frequencies in the noise. Thus, 
for sake of clarity, all the results reported in the following 
have been obtained by applying with gaussian white noise 
on Ji. 



5 Stationary solutions 

In the previous section, we built a relatively simple model 
where the cloud is described by a set of three equations 
that do not depend explicitly on the time. Such a sys- 
tem is known to be able to exhibit a complex dynamics, 
including chaos, which could explain the experimental in- 
stabilities. To know if such a complex dynamics occur in 
our conditions, the first step is to evaluate the stability 
of the stationary solutions. The stationary solutions (V s , 
Z s , N s ) of the model with the shadow effect are easily 
deduced from Eg 1 131 when the left side is put to zero. Eq. 
I13al gives immediately V s = 0. From Eq. I13cl one finds 
that the stationary solutions Z s of the CM location and 
N s of the population are linked by the simple expression: 



N s = l- zi 



(14) 



Therefore, the discussion is reduced to the analysis of Z 8 . 
Z s is the solution of Ft = in Eas ll3bl This equation can 
be resolved numerically: its global shape is illustrated in 
Fig. [21 where it is plotted as a function of Aq and no- 

A first characteristic of the Z s diagram is that Z s goes 
to 1 (i.e. N s goes to zero) at resonance. This vanishing 
of the cloud, also observed in the experiments, is a well- 
known consequence of the inefficiency of the Doppler cool- 
ing close to resonance. It is here enhanced by the shadow 
effect and the displacement of the cloud. The disappear- 
ance occurs suddenly for small no, and becomes softer as 
no increases. An interesting point is that for small no, the 
abrupt increase of Z s is linked to a very narrow bistable 
cycle. As no is increased, the bistable cycle shifts towards 
smaller Z a and smaller Aq, and the vanishing of Z s be- 
comes progressive. 
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Fig. 8. Stationary solutions of ea uat ions 1 1 3 I versus no and Aq. 
The figure represents Z a . Other parameters correspond to the 
set #1 given in tableU The definition of the N, F, SN and SF 
zones (each corresponding to different level of greys) is given 
in the text. 

The main characteristic of the Z s diagram is the pres- 
ence of several abrupt slope changes in the stationary solu- 
tions, leading to a fold in the parameter space. The shape 
of the fold depends on the parameters, in particular on 
no . Fig. E] shows four examples corresponding to situa- 
tions leading to basically different atomic dynamics. For 
no = 0.5 x 10 s (fig.m,), Z s increases smoothly with Aq (i.e. 
N s decreases slowly). The vanishing of the cloud through 
a narrow bistable cycle is not visible on the graph, as it oc- 
curs closer from resonance. As uq increases, the bistable 
cycle appears for smaller Z s (and thus larger N s ), and 
becomes physically significant. Fig. Et> shows Z s for no 
= 2.5 x 10 8 and a bistable cycle for -0.3 < A < -0.25. 
If no is further increased, the bistable cycle disappears, 
but it remains a fold corresponding to two abrupt slope 
changes of Z s versus A (fig. Et ; n — 3.4 x 10 8 ). If n is 
still increased, the fold remains, but it becomes smoother 
(Fig. Hi for no = 4 x 10 8 ). 

The dynamics of the cloud is determined by the stabil- 
ity of these stationary solutions. In particular, if no sta- 
tionary solutions are stable, complex dynamics could be 
obtained. The stability of the above stationary solutions 
is evaluated through a linear stability analysis, which as- 
sociates to each stationary solution its three eigenvalues, 
corresponding to the stability following its three eigendi- 
rections in the 3D phase space of our model. The real part 
of the eigenvalue, corresponding to a damping rate, de- 
termines the stability (stable if negative). The imaginary 
part, when different from zero, is associated to an angu- 
lar eigenfrequency, also called relaxation frequency, which 
play a main role in the dynamics. A pleasant - and simple 
- way to describe the stationary solutions is to use their 
phase space representation, where each stationary solution 
corresponds to a fixed point with its properties depend- 
ing on its eigenvalues. Let us remember that the standard 
terminology distinguishes the stable node (all eigenvalues 
real and negative), the stable focus (all real parts negative, 
two eigenvalues complex conjugate), the saddle node (all 







(d) ; ' / 




— — ~~ _lb)___— - 




(a) 



-0.5 -0.4 -0.3 -0.2 

A 



Fig. 9. Evolution as a function of the detuning of the sta- 
tionary solution Z s of equations 1131 The full (resp. dashed) 
line corresponds to a stable (resp. unstable) solution. In (a), 
n = 0.5 x 10 8 ; in (b) n = 2.5 x 10 s ; in (c), n = 3.4 x 10 8 ; in 
(d), no = 4 x 10 s Other parameters correspond to the set #1 
of Table H 

eigenvalues real, at least one positive) and the saddle focus 
(at least one real part positive, two eigenvalues complex 
conjugate). For sake of simplicity, this terminology will be 
used in the following. 

The set of points of vertical tangency in Z s determines 
a line which delimits the unstable stationary solution of 
the medium branch of the bistability cycle. Linear stabil- 
ity analysis shows that the fixed point in that area is a 
saddle-node with real eigenvalues, two being positive and 
one negative. But the unstable stationary solutions extend 
beyond this area. In particular, the stationary solutions 
can also be unstable on a part of the upper branch of the 
bistable cycle, and even outside the bistable cycle, when 
the stationary solution is unique. This is illustrated on fig. 
|H1 where the nature of the stationary solution is indicated 
by a level of gray and a code. In the SN zone, the fixed 
point associated with the stationary solution is a saddle 
node, i.e. the stationary solution is unstable, and its three 
eigenvalues are real, two being positive. In the SF zone, 
the fixed point is a saddle focus, i.e. the stationary solu- 
tion is unstable, and two of its eigenvalues are complex 
with a positive real part, and the third one is real nega- 
tive. In the F zone, the fixed point is a stable focus, i.e. 
the stationary solution is stable, and two of its eigenval- 
ues are complex with a negative real part, the third one 
is real negative. Finally, in the N zone, the fixed point is 
a stable node, i.e. the stationary solution is stable and all 
eigenvalues are real negative. The dashed line indicates the 
location of the bifurcation, i.e. the transition from stable 
to unstable stationary solutions. In most cases, it occurs 
from stable focus to saddle focus, through a super-critical 
Hopf bifurcation. 

As the detuning is varied, four typical situations may 
occur, already illustrated in Fig. El For small n (Fig.EK), 
Z s is always stable, and its dependence versus Aq is al- 
most flat (except very close from resonance): we expect 
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a stationary cloud slightly moving with the detuning. In 
the bistability area (Fig. Eb), the central branch of the 
bistability cycle is unstable, as usual in such a situation, 
but the upper branch is also partly unstable, reducing 
drastically the parameter range where the system is ef- 
fectively bistable. However, a narrow bistable zone should 
be observed. As no is still increased, only one stationary 
solution remains (Figs Eb): which is unstable on the fold. 

Finally, for large no (Fig. [IJO, the solution is always 
stable (F zone), except very close from resonance. Thus 
we expect to observe a stationary cloud moving with Aq. 
However, because of the fold, this moving is non linear, 
and on the fold, for —0.45 < A < —0.35, Z$ is more 
sensitive to the Aq value, because the slope here is larger. 
On the other hand, a closer look to the eigenvalues (Fig. 
I10|) shows that at the level of the fold, the real part A of 
the complex eigenvalues A ± iuu is close to zero, while the 
eigenfrcquency crosses a minimum. This is more evident 
on Fig. ^2 where the same situation is illustrated for a 
different set of parameters. The vanishing damping rate 
|A| associated to an increased sensitivity to the parameter 
values means that a small perturbation will cause a large 
undamped reaction of the system, with relaxation oscil- 
lations at a frequency oj crossing a minimum while Aq 
is varied on the fold. In the fold inflection point, where 
the slope is maximum, the eigenfrequency and damping 
rate are minimum: the effects of perturbations is expected 
to be maximum at this point. Note that the evolution of 
A and oj around the inflection point is asymmetric: the 
damping rate and relaxation frequency increase rapidly 
for detunings smaller than the inflection point, while they 
remain of the same order of magnitude on the other side. 
This asymmetry is expected to have consequences on the 
dynamics. 




Fig. 10. Evolution as a function of the detuning of the station- 
ary solution Zs and its eigenvalues, for the parameters of Fig. 

The stationary solution is given through the full bold line 
in the upper side of the figure. The dashed line noted oj repre- 
sents the imaginary part of the complex eigenvalues, while the 
full line corresponds to their real part A. The third eigenvalue 
is always real negative, almost constant and everywhere larger 
than A. It cannot be distinguished from the zero axis at the 
scale of the figure. 

Deterministic instabilities are expected to occur for pa- 
rameters where all stationary solutions are unstable, i.e. 
in the SF monostable zone^Dl- In ai l other areas, the sta- 
tionary solutions are stable, and therefore, deterministic 




Fig. 11. Same as Fig.^Jfor parameter set #2 of Tabled For 

sake of clarity, the evolution of A has been truncated for small 
detunings, where it becomes almost linear, reaching e.g. the 
value of A = -1500 for A = -0.3. 

instabilities cannot occur. However, the presence of the 
fold and the proximity of a Hopf bifurcation generate a 
type of stochastic behavior similar to instabilities^- 

Indeed, on the one hand, the proximity of a Hopf bi- 
furcation is known to favor the appearance of coherence 
resonance JOi., as discussed in j^]. Let us recall that coher- 
ence resonance is the counterpart of stochastic resonance 
in autonomous systems |1 5). Stochastic resonance appears 
in some forced systems: it may be seen as an amplifica- 
tion by noise of the system response to a forcing. In other 
terms, the signal to noise ratio of the output periodic sig- 
nal resulting from the modulation exhibits a maximum 
when the noise amplitude increases. In stable autonomous 
systems, the behavior is not periodic, but a phenomenon 
similar to stochastic resonance can lead to an amplifica- 
tion of an internal resonance: it is the internal stochastic 
resonance, or coherent resonance (1 71118] . Coherent reso- 
nance is also known to generate noise induced coherent 
oscillations, which explain some of our experimental ob- 
servations, as shown in !) . 

On the other hand, as shown above, the particular con- 
figuration of the eigenvalues on the stable fold makes the 
system very sensitive to small perturbations, i.e. to noise, 
and thus we expect a large increasing of the response to 
noise on the fold. The resulting behavior looks like in- 
stabilities, whereas it is simply noise amplification. The 
behavior obtained on the stable fold when noise is taken 
into account is detailed in the next section. 



6 The stable fold: stochastic instabilities 

In ^3 > it has been shown that the model is able to repro- 
duce the C instabilities, in the monostable SF zone. On 
the contrary, S instabilities do not appear in this situa- 
tion. In fact, it has been shown in jH] that S instabilities 
are not instabilities in the usual significance, but the re- 
sult of the amplification of the system intrinsic noise, and 
thus it is necessary to add noise in the simulations to ob- 
serve S instabilities. This dynamics could appear when the 
stationary solution of the MOT is not stable, but would 
be difficult to analyze, because in this case the resulting 
dynamics would be the superimposition of the determinis- 
tic instability with the stochastic motion, the latter being 
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masked by the former. For this reason, the present sec- 
tion is devoted to the study of the influence of noise on 
the stable stationary solutions of the model, in particular 
in the vicinity of the fold. As the model here is slightly 
different from that used in HJ, we could expect different 
results from those obtained in It is shown in the fol- 
lowing that the modifications introduced in the model do 
not alter the previous conclusions, i.e. that S instabilities 
are produced by noise amplification on the fold. 

To evaluate the response of the system to noise in the 
vicinity of the fold, we have plotted the amplitude of the 
motion perturbed by noise, versus the detuning, across the 
fold (Fig. 112(1 . The amplitude of the motion is measured 
through the standard deviation a of Z. Fig. ^] has been 
obtained for a gaussian white noise applied on I\ with an 
amplitude of 10 -3 , in a situation where the stationary so- 
lutions are stable on the fold. Noise amplification appears 
clearly on the fold, and mimics instabilities, as in the ex- 
periments (Fig. The maximum of the motion ampli- 
tude, obtained for A = —0.23, corresponds to a standard 
deviation of 1.1 x 10~ 3 , i.e. 33 /im, in good agreement 
with the experimental value. By comparison, the motion 
amplitude is 0.036 for Aq — —3, i.e. the effect of noise is 
30 times larger on the fold than in Aq = —3. Note that, as 
similar motion amplitudes are obtained in the experiment 
and in the model for a noise level £ = 10~ 3 , this gives us 
in return an evaluation of the experimental noise level. 




-0.30 -0.28 -0.26 -0.24 -0.22 -0.20 -0.18 

A 



Fig. 12. Evolution of the signal standard deviation a (full line) 
and of the cloud location stationary solution Zs (dashed line) 
as a function of the detuning Aq. White noise is applied to 7i 
with £ = 10 -3 . Other parameters correspond to the set #2 of 
Tab.il 

Thus, as for the model developed in [S], the fold plays 
the role of a noise amplifier, leading to an unstable cloud 
over a limited range of the parameters. The maximum of 
the motion amplitude appears at the fold inflection point. 
The amplitude decreases progressively and quasi symmet- 
rically on each side of this point, in spite of the asymme- 
try of the eigenvalues discussed in the previous section. 
The range where the instabilities can be observed, ap- 
pears to be of the order of the "width" of the fold, that 
we could define as the interval between the two abrupt 
slope changes delimiting the fold (a more precise defini- 
tion could be given, but is not useful for the following). 
In the case discussed here, the range in detuning is 0.02, 



i.e. much smaller than the experimental range, which was 
typically of the order of 1. However, the width of the slope 
is very dependent on the MOT parameters: for example, 
on Fig.E3 where I\ and no are different, the fold is more 
than twice larger than in the present case. Thus there is 
no doubt that it is possible to find a parameter set giv- 
ing a correct fold width. However, because of the extreme 
simplicity of the model as compared to the experiments, 
such a search has no meaning: the aim here is just to show 
that very few ingredients, including the shadow effect and 
the noise, are able to reproduce the global experimental 
behavior. We have now demonstrated that the noise re- 
sponse can effectively have the appearance of instabilities 
for some parameters. We will now examine in detail the 
resulting dynamics, and compare it to the experimental 
observations. 

Fig. shows the time evolution of Z and N in Aq = 
—0.23, where the amplification is maximum. The signal 
is of course stochastic, but several differences appear, on 
the one hand between the cloud response and the applied 
noise, and on the other hand between the Z and N evo- 
lutions. Indeed, while the applied noise is white, a fre- 
quency dominates in the response, in particular for Z. 
This is clearly visible on the signal, and of course on the 
spectrum (Fig. 114(1 , where a peak appears at a frequency 
uis — 38 Hz. For N, while the frequency lus still appears, 
low frequency fluctuations dominates, as it is also shown 
by the spectrum fFig. 114(1 . This behavior with two differ- 
ent characteristic times is very similar to the experimental 
Sh behavior illustrated in the Fig. H2 in both cases, the 
higher frequency component dominates the motion of the 
atoms, while the number of atoms in the cloud is mainly 
driven by the low frequency component. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

time (s) 



Fig. 13. Dynamics obtained by resolving Eqs. 1131 . with the 
parameter set #2 and Ao = —0.926. White noise is applied to 
ii with C = 10" 3 . 
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Fig. 14. Power spectra obtained by resolving Eqs. ll3l with the 
same parameters as in Fig. 1131 



In the experiments, the origin of the higher frequency 
was not identified, as it does not correspond to a known 
experimental characteristic frequency. However, we know 
now that the stable stationary solution on the fold has 
complex eigenvalues, and thus is associated with a relax- 
ation frequency. As noise is known to be able to excite such 
non linear resonance eigenfreauencies|19|. it is interesting 
to compare in the simulations the frequency appearing in 
the dynamics with the eigenfrequency of the stationary 
solution. This is illustrated on Fig. EJf° r different values 
of the detuning. The full and dashed lines concern the sta- 
tionary solutions: they shows the evolution of ui and Z$ 
versus A . They are a reproduction of Fig. El an d are re- 
called just for comparison. The squares give the resonance 
frequencies of the noisy dynamics. They are obtained by 
fitting the calculated spectra of Z, as illustrated on Fig. 
114b . with a lorentzian function. Similar fits on the TV spec- 
tra give the same results. Fig.ll5lshows that for detunings 
larger than the fold inflection point, the noisy resonance 
frequency and the eigenfrequency correspond exactly. This 
confirms that the frequency appearing in the dynamics is 
the relaxation frequency appearing in the eigenvalues as- 
sociated with the stationary solution. Thus "instabilities" 
appear in fact as a noisy amplification of the relaxation 
frequency of the MOT, through a phenomenon already 
observed in many other systems, as e.g. in lasers the 
small damping rate allows noise to excite the relaxation 
frequency, altering the frequency distribution of the sys- 
tem response to noise. 

For detunings smaller than the fold inflection point, 
Fig. El exhibits differences between the eigenfrequencies 
and the resonance frequencies. A close look at the dynam- 
ics of the cloud shows that the main point is not this dif- 
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Fig. 15. Evolution of the resonance frequencies as a function 
of the detuning. The full and dashed line recalls the evolution 
of ui and Zg, already shown on Fig. 1111 The squares represent 
the resonance frequencies appearing in the dynamics, obtained 
by fitting a lorentzian function on spectra similar to those of 
fig. [H The bars give the width of the lorentzian. 



ference, but a dramatic broadening of the noise resonance, 
leading to an indetermination of the resonant frequency. 
This is illustrated on Fig. 1151 where the bars associated 
with each point have a height of Av, the width of the 
lorentzian peak obtained by fitting the Z spectrum. On 
the inflection point, the resonance is narrow, as it appears 
on Fig. El O n the right of the inflection point, the reso- 
nance remains narrow, reaching a maximum of Av = 14 
Hz in Aq = —0.2. On the contrary, for detunings smaller 
than the inflection point, Av increases rapidly, reaching 
already the value Av = 30 Hz in Aq = —0.235 (the in- 
flection point is in Aq = —0.23). It is clear on Fig. El 
that the discrepancy between the noise resonance and the 
relaxation frequencies is linked to the width of the reso- 
nance: as the resonance broadens, it also flattens, and the 
central frequency becomes irrelevant. This is the reason 
why the resonant frequencies for detunings smaller than 
Aq = —0.26 are not reported on the figure. 

This asymmetry in the dynamics around the fold in- 
flection point is of course linked to the asymmetry in the 
eigenvalues associated with the stationary solution. For 
detunings larger than the inflection point, the damping 
rate remains small (less than 45 s _1 ), and thus the noise 
excitation of the relaxation oscillations remains efficient. 
On the contrary, for detunings smaller than the inflec- 
tion point, the damping rate increases rapidly, decreasing 
the efficiency of the excitation, and leading to a flat reso- 
nance. The consequence on the dynamics is illustrated on 
Fig. El where the power spectra of the Z dynamics are 
represented for two different values of the detuning, lo- 
cated symmetrically with respect to the inflection point. 
In (a), for Aq = —0.22, i.e. for a larger detuning than 
the inflection point, the resonance, centered on 50 Hz, re- 
mains narrow (Av = 11 Hz), and thus significant. In (b), 
for Aq = —0.24, i.e. on the other side of the inflection 
point, the resonance, centered on 100 Hz, is already 60 Hz 
wide, and also five times lower (note the different vertical 
scales in (a) and (b)). Such a resonance is no more sig- 
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nificant from an experimental point of view: indeed, the 
experimental dynamics is analyzed from time series with 
a necessary limited number of points, leading to spectra 
with a resolution much smaller than in the simulations. It 
is clear that a resonance as those observed on the left of the 
inflection point could not be detected in the experiments, 
and thus one expects that in this situation, experiments 
deliver non resolved spectra. It was effectively the case for 
Sl instabilities (Fig.0J), and thus we can conclude that the 
Sl dynamics observed in the experiments is also explained 
by the present model. 
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Fig. 16. Z spectra of the instabilities for two values of the 
detuning situated symmetricallty with respect to the fold max- 
imum. In (a), A = -0.22; in (b), A = -0.24. Note that the 
vertical scale in (a) is five times larger than in (b). 

The results obtained in this section can be summa- 
rized in the following way: when the parameters of the 
cloud are such that the unique stationary solution is sta- 
ble, and therefore that no deterministic instabilities can 
occur, the action of noise mimics instabilities. This origi- 
nates in a fold of the stationary solution, which plays the 
role of a noise amplifier. Moreover, the particular prop- 
erties of the eigenvalues associated with the stationary 
solutions on the fold - small damping rate and relaxation 
frequency - lead to the existence of the stochastic regime: 
for detuning larger than the fold inflection point, a reso- 
nance is excited by the noise, leading to a dynamics with 
a dominating frequency; for detuning smaller than the in- 
flection point, the damping rate increases, and the reso- 
nance vanishes. This allows us to interpret all the regimes 
observed experimentally. 



The particular properties of the eigenvalues discussed 
above are in fact linked to the proximity of a Hopf bifur- 
cation. Another known consequence of such a situation is 
to favor the existence of coherence resonance, and such a 
phenomenon exists effectively in our case Now, an in- 
teresting question is to determine what are the exact con- 
nections between all these phenomena, and in particular 
if the existence of coherence resonance has a determining 
role in that of instabilities. The answer is obviously: no! 
Indeed, the essential ingredient for noise amplification is 
the fold; the properties of the eigenvalues only determine 
the time characteristics of the dynamics. In particular, we 
have checked that with the present model, no coherence 
resonance occurs when the MOT exhibits Sl instabilities. 

7 Conclusion 

It has been shown recently that the cloud of cold atoms 
obtained from a magneto-optical trap may exhibit a com- 
plex dynamics in the regime of high atomic densities. The 
observed behaviors can be essentially separated in two dif- 
ferent types, depending of their nature: stochastic insta- 
bilities 9 or deterministic instabilities 10 . The aim of this 
paper was to describe extensively the experimentally ob- 
served stochastic dynamics, and to understand its mech- 
anisms, through a model briefly presented in [10) and de- 
tailed here. 

A detailed analysis of the experimental results in the 
stochastic regime shows a variety of dynamical behaviors, 
which differ by the frequency components appearing in 
the dynamics. Indeed, some instabilities exhibit only low 
frequency components, while in other second time 

scale, corresponding to a higher frequency, appears in the 
motion of the center of mass of the cloud. 

The simple stochastic ID-model that we use here al- 
lows us to retrieve and interpret these experimental dy- 
namics. The model shows that the existence of instabil- 
ities is linked to folded stationary solutions where noise 
response is enhanced. Moreover, the proximity of a Hopf 
bifurcation and the resulting conditions on the stability 
of the stationary solutions - small damping rate and ex- 
istence of a relaxation frequency - explains the existence 
of several types of regimes: indeed, depending on the pa- 
rameters, noise is sometimes able to excite the relaxation 
frequency, leading to the appearance of the second time 
scale in the dynamics. Globally, the agreement between 
this over-simplified ID model and the 3D experiments is 
surprisingly good. Not only it allows us to make a qualita- 
tive interpretation of the experimentally observed dynam- 
ics, but also gives quantitative results compatible with the 
experimental measures, in particular concerning the am- 
plitude of the instabilities. However, it is clear that a 3D 
model is necessary for an accurate quantitative description 
of all the experimental observations. 

The model also emphasized the close relations between 
the stochastic and deterministic instabilities. Indeed, both 
types of dynamics appear to be link to the same factors: 
the effective dynamics depends mainly on the distance be- 
tween the working point and the bistable cycles (or the 
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Hopf bifurcation). In fact, the dynamics described in the 
present paper appears as the first sign of the deterministic 
instabilities described in ^U], and thus may be considered 
as a noisy precursor to deterministic instabilities 8 , 20 . A 
study of the phenomenon from this point of view should 
put in evidence new properties of these instabilities. 
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